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ABSTRACT 

Star clusters are known to have smaller intrinsic metallicity spreads than 
dwarf galaxies due to their shorter star formation timescales. Here we use in- 
dividual spectroscopic [Fe/H] measurements of stars in 19 Local Group dwarf 
galaxies, 13 Galactic open clusters, and 49 globular clusters to show that star 
cluster and dwarf galaxy linear metallicity distributions are binomial in form, 
with all objects showing strong correlations between their mean linear metallic- 
ity Z and intrinsic spread in metallicity o(Z)'^ . A plot of cr(Z)^ versus Z shows 
that the correlated relationships are offset for the dwarf galaxies from the star 
clusters. The common binomial nature of these linear metallicity distributions 
can be explained with a simple inhomogeneous chemical evolution model (e.g.. 



Oey (2000)), where the star cluster and dwarf galaxy behaviour in the cr(Z)^ — Z 
diagram is reproduced in terms of the number of enrichment events, covering 
fraction, and intrinsic size of the enriched regions. The inhomogeneity of the 
self-enrichment sets the slope for the observed dwarf galaxy cr(Z)^ — Z corre- 
lation. The offset of the star cluster sequence from that of the dwarf galaxies 
is due to pre-enrichment, and the slope of the star cluster sequence represents 
the remnant signature of the self-enriched history of their host galaxies. The 
offset can be used to separate star clusters from dwarf galaxies without a priori 
knowledge of their luminosity or dynamical mass. The application of the inhomo- 
geneous model to the oiZ^' — Z relationship provides a numerical formalism to 
connect the self-enrichment and pre-enrichment between star clusters and dwarf 
galaxies using physically motivated chemical enrichment parameters. Therefore 
we suggest that the oiZ^' — Z relationship can provide insight into what drives 
the efficiency of star formation and chemical evolution in galaxies, and is an 
important prediction for galaxy simulation models to reproduce. 
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Subject headings: galaxies: dwarf - galaxies: abundances - globular clusters: general 
- open clusters and associations: general 
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Introduction 



Correlations between the average metallicity and the mass or luminosity of a galaxy 



e.g., Lequeux et al. 1979 Skillman et al. 1989), the mass-metallicity relation (MZR), offers 



a window into how chemical evolution in galaxies proceeded over cosmic times. Numerous 
studies have explored the interplay between chemical enrichment produced via SNe versus 



the efficiency of outflows - and in recent years gas inflows (Keres et al. 2005 Brooks 



et al. 2007, 2009). However, a complete picture for how the MZR came to be at redshift 



zero, and evolved before that, is still missing at this time. The star formation (SF) and 
metal retention efficiencies as a function of galaxy mass are also needed to understand the 
distribution of baryons in dark matter halos. For example, to what degree is SF modulated 



by metallicity dependent H2 formation (Kuhlen et al. 2011), heating from stellar and SN 



feedback (Governato et al. 2010), or environment (Mayer et al. 2006). It is also not clear 
why galaxies with masses less than ~ 10^'^Mq have such low baryon fractions relative to 



Milky Way (MW) sized galaxies (e.g., Guo et al. 2010) 



Along with the mean metallicity of a host system, comparisons of the intrinsic spread 
in metallicity for different objects may help constrain the timescale for chemical enrichment 
in galaxies. There is increasing observational evidence for inhomogeneous mixing in dwarf 



galaxies (e.g., Koch et al. 2008 Venn et al. 2012), or alternatively hierarchical merging of 



minihalos (Salvadori & Ferrara 2009, 2012), both of which can be traced by the dispersion 



in elemental abundances. 



Kirby et al. (2011) and Norris et al. (2010) each looked at the spread in [Fe/H] versus 



host galaxy luminosity (L) amongst Local Group dwarf spheroidals (dSphs) and ultra faint 
dwarfs (UFDs), with both studies finding a clear anticorrelation. Leaman et al. (2012b) 
included higher luminosity dwarf irregular (dirr) systems and found that the spread in 



[Fe/H] tended to plateau above L = IO^Lq. Alternatively, Kirby et al. (2011) showed that a 
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strong correlation between spread in linear metal fraction, cr(^) and luminosity exists, with 
much tighter scatter than in the a{[Fe/ H]) — L plane. It is currently unclear whether these 
trends in metallicity spread are due to SF duration or enrichment efficiencies. Nevertheless 
such correlations may be useful to infer the relative formation and enrichment timescales 



for the UFDs and dSphs (i.e., Willman et al. 2011). In fact, the similarities in metallicity 
spreads between the UFDs and the more luminous dSphs is one item used to differentiate 



the UFDs from simple star clusters (Willman & Strader 2012) 



Both star clusters and dwarf galaxies (DGs) have resolved spectroscopic observations 
of many stars, however they are treated independently in chemical evolution studies. 
Specialized n-body+SPH simulations exist for both dwarf galaxies and star clusters (e.g.. 



Governato et al. 2010 Nakasato et al. 2000), however they are not general enough, with 



sufficiently high resolution, to analyze both objects in a single simulation (although see 



Griffen et al. 2010). Therefore, comparative chemical analyses of star clusters and dwarf 



galaxies have been typically done through analytic descriptions. Chemical evolution models 
that relax some of the homogeneity and mixing assumptions from the simple one-zone 
closed box models have also been developed, and these are able to parameterize physical 



processes (such as stellar feedback and outflow) in low mass systems (e.g., van den Hoek & 



de Jong 1997 Argast et al. 2000 Karlsson 2005). These chemical evolution models also 



incorporate stochastic sampling effects, allowing for detailed predictions in the elemental 
abundance distributions. 

Specialized inhomogeneous mixing models can generate simulated metallicity 
distribution functions (MDFs), however comparisons of these to observed MDFs are very 
sensitive to biases in the sampling, radial coverage, and measurement methods. In this 
paper, we use an inhomogeneous stochastic chemical evolution model to compare only the 
intrinsic spread and mean metallicity to observationally well sampled star clusters and 
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dwarf galaxies. This analysis does not require an estimate of the host galaxy luminosity (a 
quantity which may be strongly affected by tidal stripping), and exploits specific statistical 
properties of the MDFs which we demonstrate are universal. The scaling relations developed 
here show the link between the formation of star clusters and chemical evolution of dwarf 
galaxies, through a simple binomial distribution in metal enrichment events. 



2. Data Sources 

The data for this paper come from published spectroscopic studies of individual stars 
in Local Group dwarf galaxies, open clusters, and globular clusters in the MW. [Fe/H] 
values and errors were used to compute metallicity distributions and intrinsic dispersions 
uniformly. The metallicity distributions of the star clusters are known to exhibit small 



spreads in [Fe/H] (Carretta et al. 2009) and were initially included only to offer an 



interesting benchmark to compare to the dwarf galaxy distributions and dispersions. 

In assembling the data, an effort was made to use homogeneous studies of several 
objects observed with a single instrument setup and analysis technique. However in cases 
where a choice of studies were available, data were selected where: 1) spectral resolving 
power was the greatest, 2) the number of stars in the study was large, 3) in the case of 
Calcium triplet (CaT) based [Fe/H] measurements, the equivalent width data of individual 
stars was available allowing us to transform these to [Fe/H] using the improved non-linear 



calibration by Starkenburg et al. (2010), and 4) the spatial completeness was high. For 
some systems, such as the UFDs, the small number of stars is due to the intrinsic low 
luminosity of the system. These small number statistics must be kept in mind as the 
sample distributions may deviate from the true system distributions, especially in the case 
of spatial incompleteness and radial population gradients (c.f., Leaman et al., 2012b). We 
restrict ourselves to spectroscopic abundances to avoid the age-metallicity degeneracies 
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which plague photometric metallicities. Table 1 lists the objects considered in this Paper, 
along with the original literature study, number of stars with spectroscopic metallicity 
measurements, resolving power of the spectrograph used, and [Fe/H] measuring method - 
CaT, or medium resolution/high resolution spectroscopy (MRS/HRS). 

Spectroscopic measurements of individual stars are used to derive [Fe/H] abundances by 
comparing synthetic spectra to observed Fe lines, by measuring the equivalent width of Fe I 
and Fell lines, or relating the Calcium Triplet equivalent width to [Fe/H] using empirical 
relations (i.e., Armandroff & Da Costa 1991 Starkenburg et al. [2010 ). The individual 
[Fe/H] measurements are taken directly as listed in each of the literature studies, with a 
few exceptions: In the case of high resolution abundances, the [Fe/H] measurements were 
taken as the weighted average of the [Fel/H] and [Fell/H] values, however their variation 
was typically less than 5%. In surveys using the CaT, [Fe/H] values were recomputed using 



the two strongest CaT lines and the non-linear relations of Starkenburg et al. (2010). This 



recalibration ensures that the non-linear bias at faint magnitudes and low metallicities 



is removed. Another exception is the open cluster survey by Warren & Cole (2009) and 



Cole et al. (2004) where the V or K band magnitude calibrations (created in those papers) 
were used. The near solar metallicity of the open clusters means that [Fe/H] biases are 
negligible, unlike in the lower metallicity globular clusters or dwarf galaxies. Finally, the 
high resolution studies may have used different solar abundances for their zero points (e.g., 
logeFe ~ '''•52 from ?, or 7.45 from ?), however the differences are very small, having no 
significant effect on the results presented in this paper. 



3. Metallicity Distributions 



With a relatively homogeneous set of [Fe/H] values for the dwarf galaxies, open 
clusters, and globular clusters, metallicity distribution functions can be computed for all 
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objects. This was first done in [Fe/H] space, and subsequently, each individual stellar [Fe/H] 
measurement was transformed into a linear metal fraction Z, assuming; 
This conversion has two sources of uncertainty; the first due to the exclusion of contributions 
from a-elements to the total metal fraction. At low metallicities the a element contribution 



is larger, but is roughly constant with metallicity (e.g., Venn et al. 2004). Variations in 
[a/Fe] are noticeable (< 0.4 dex) at intermediate metallicities, particularly in dwarf galaxies 
(e.g., Tolstoy et al. 2009), however inclusion of [a/Fe] in each galaxy is difficult and beyond 
the scope of this paper. Therefore, while a true estimate of the total linear metal fraction 
is not encompassed in the Z value, we use Z to track the contribution of Fe as a simple 
"linear [Fe/H]" or "linear metallicity" in order to study the metal dispersions in a linear 
scale. The second source of uncertainty is in the solar value of Zq, which may be as high 



as 50% (VandenBerg et al. 2007). The adopted Zq values from the literature represents a 



small offset that may contribute to the scatter in the relations discussed in §4. 

The metallicity distribution functions in [Fe/H] and Z are shown in the left and middle 
panels of Figure [l] for four representative dwarf galaxies. While the [Fe/H] distributions 
are relatively Gaussian, the linear metallicities follow a Poisson, or binomial, distribution 
(especially noticeable in the lower metallicity systems such as Sextans). Comparing the 
linear metallicity spread (j{Z) to metallicity Z or other galaxy properties must be done with 
care when the distributions are Poisson or binomial, as the variance of such distributions 
is proportional to the mean metallicity of the object, and physical comparisons between 
objects may be difficult to disentangle from the statistical properties of the distribution 
itself. Thus, it should be kept in mind that spatial biases in a sample will influence the 
spread in [Fe/H] and Z differently. For example, a galaxy with a steep negative metallicity 
gradient which is sampled preferentially in the inner regions will return a sample mean 
which is obviously higher than the total galaxy for both [Fe/H] and Z, however the intrinsic 
spread in Z will be decreased less than the spread in [Fe/H] due to the tail of the Z 
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distributions being in the metal-rich regime. 

3.1. Quantifying the Statistical Form of the Z Distributions 

There is a strong correlation between mean metallicity and intrinsic metallicity spread 
among the objects, as seen in the right hand panels of Figure [1} This correlation can be 
quantified by looking at the relative dispersion index (defined as D = cr(Z)^/Z) for each 
of the galaxies and star clusters. This index can be used to infer the level of correlation 
between the variance and mean of a Z distribution; for a Poisson distribution, D = 1, while 
distributions with D > \ and D < \ are described by a negative binomial and binomial 
distribution, respectively. The dwarf galaxies and star clusters all show dispersion indices 
less than 1, and are therefore described by binomial distributions. This observational result 
is interesting because these objects share the same functional form in their linear metal 
distributions, despite the range in masses and environments. 

The globular clusters typically show the lowest values of ranging from 10^^ < 
D < IQ-^, while the more luminous dwarf galaxies show distributions which are closer 
to a Poisson distribution, with the dispersion indices ranging from 10"^ < D < 10"\ 
Generally a binomial distribution can be thought of as a set of n trials that occur with 
some probability of success q (see §5.1 below). This is a specific case of the more general 
Poisson-binomial distribution, which is a sequence of n trials which each have an individual 
probability of success of g^. In the limit that all the probabilities are equal, then the 
Poisson-binomial reduces to the standard binomial distribution. The interesting question of 
why the distributions of linear metallicity all exhibit binomial forms will be discussed in §4. 
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4. Intrinsic Metallicity Spreads 

To analyze the intrinsic spread in metallicity for the objects in Table 1, we begin 
by removing the contribution to the metallicity spread due to observational errors. To 
calculate the intrinsic metallicity spreads (j{[Fe/H]) and <j{Z), we follow the description by 



Kirby et al. (2011) to solve the following equation numerically for a{[Fe/ H]): 

N 



N ^ 



i[Fe/H], - {[Fe/H])f _ ^ 



N^^{6[Fe/H],y + a{[Fe/H]y 

Where [Fe/H]i are the individual spectroscopic measurements from each star, and S[Fe/H]i 
are the individual measurement errors per star. An analogous equation is used to solve 
for the intrinsic linear metallicity spread cr(Z), with the the error in linear metal fraction 
computed from the [Fe/H] measurements and errors as: 

^ = (Z,/Zo)ln(10)5[Fe/H]„ (2) 

Errors on the intrinsic spreads are computed by numerical jackknife estimates, where 
each star is removed from the sample in turn, and the mean and variance recomputed. The 
error on the mean or variance is then, 



N -1 

where y is the quantity (mean, or variance) of interest. In all cases, an upper limit to 
the intrinsic metallicity spread o"([Fe/H]) has been assigned to the distribution when the 
uncorrected metallicity spread has a smaller value than the smallest measurement error of 
an individual star, i.e., the upper limit value is taken as the smallest measurement error. 
Table 2 shows the computed average metallicities and intrinsic metallicity spreads for the 
objects we consider. 
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4.1. Consequences of the Binomial Nature of the Z Distributions for 
Interpreting Intrinsic Metallicity Spreads 

To visualize the binomial nature of the linear Z distributions in the dwarf galaxies and 
clusters, Figure [2] plots the intrinsic variance, ct(Z')^, versus the mean linear metallicity, 
Z. There is a correlation between the mean metallicity and intrinsic spread in metallicity 
for both the dwarf galaxies and star clusters, characteristic of a binomial distribution as 
discussed in §3.1. A clear separation in the sequences is evident, as dwarf galaxies at 
all metallicities show intrinsic metallicity spreads larger than the clusters at a given Z . 
The RMS scatter about the dwarf galaxy (DG) line is only 0.10 dex, while it is more 
than 7 times larger about the cluster sequence (0.71 dex). Interestingly however, the star 
clusters do not scatter up and overlap with the galaxy population at any point. This 
discrete offset is related to the fact that the star clusters have been pre-enriched within 
systems which themselves lie on the galaxy sequence. Therefore the slope of the star cluster 
sequence, having the observational errors accounted for, represents a remnant signature of 
the self-enriched history which their formation environment (galaxies) underwent (see §5.3). 

The split in the DG and star cluster sequences is robust to both sampling errors and 
internal biases, i.e., if the quoted errors (5[Fe/if]i) in the original studies are under or 
overestimated. To illustrate this point, in the right panel of Figure [2] the cr(Z)^ — Z diagram 
is shown with the error on each object from the jackknife tests. Objects with very few 
stars in them, such as Willman 1 and Segue 1, show larger errors, but still remain offset 
from the star cluster sequenc^ The open circles in Figure |2| show the change in the 

^We stress that this test is not sufficient to prove that Willman 1 and Segue 1 are dwarf 
galaxies since one or two Milky Way foreground stars in these small systems would enlarge 
their metallicity spreads so that they lie on the DG sequence. Proper foreground removal 
using statistical photometric and kinematic techniques is still crucial to understanding the 
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intrinsic dispersion if the reported uncertainties on the [Fe/H] values, (5[Fe/H], were under 
or over estimated by 50 — 150% in 10% increments. The sequences still show no overlap 
and there are no apparent correlations between the number of observed stars or quality of 
the observations. The cr(Z)^ — Z appears to be a robust way to differentiate star clusters 
from dwarf galaxies, and therefore we suggest that this as a new and direct method for 
distinguishing these stellar populations. 

It is important to note that the correlation of intrinsic linear metallicity spread oiZ)^ 
with average metallicity Z is more nuanced than just assuming this is the natural result of 
how the metallicity increases with age of a system. While systems with longer epochs of 
star formation will have larger spreads in metallicity, the DG sequence in Figure [2] is not a 
sequence in SFH (e.g., the positioning of Willman I is close to Fornax). Similarly, if the 
expression of SFH were the dominant factor determining the correlations in Figure |2] then 
one might expect such a correlation between metallicity spread and mean metallicity to 
hold in the [Fe/H] measurements. However, as we show in Figure [sj the intrinsic spread 
a{\Fel H\) is independent of the average [Fe/H] of the dwarf galaxies. Thus the usefulness 
of the oiZ^' — Z diagram can be seen - this relationship is due to the discrete binomial 
form of the chemical evolution and linear metallicity distributions in these systems, and it 
is not simply a product of generic metal enrichment. 

4.2. Revisiting Metrics for Separating Star Clusters and Dwarf Galaxies 

With the intrinsic value of the metallicity spread, a{\Fej H\) or and the associated 
errors, we compare the behaviour of the metallicity dispersions of the star clusters, and 



nature of these systems, e.g., Aden et aL 2009 McConnachie fc Cote 2010 [Koposov et al. 
20TTI IMunoz et al. 1120 12 
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dwarf galaxies in the a{[Fe/H]) — L and cr{Z) — L planes in Figure 111 As expected, the 



clusters show very low spreads in metallicity relative to the dwarf galaxies (c.f., [CarrettaT 



et al. 2009 Willman & Strader 2012). However in the top panel of Figure 4 the high 



luminosity systems show a fiat trend in a{[Fe/H]) — L, and as expected the M54/Sgi[^and 
uCen are the only potential clusters that show a large spread in heavy metals - consistent 
with the theory that both may be the nuclear remnants of accreted dwarf galaxies (e.g., 



Bekki & Freeman 2003 Carretta et al. 2010). 



In the bottom panel of Figure |4] the dwarf galaxies show a tight correlation in (j{Z) — L, 
as first reported by Kirby et al. ( |2011[ ). Interestingly, the globular clusters and open clusters 
do not show this trend, and yet the former overlap in the parameter space of the dwarf 
galaxies with luminosities of 10^ — IO^Lq. The linear metallicity spread of an object in this 
luminosity range is then not necessarily indicative of it being a galaxy. 

Given the well known luminosity-metallicty relation exhibited by Local Group dwarfs 



(Dekel & Woo 2003 Kirby et al. 2011), plotting a{Z) versus L is equivalent to o-{Z) versus 



average metallicity Z. As discussed in §3 the Z distributions are clearly non-Gaussian, and 
in fact they are binomial. Therefore it is not surprising that a tight sequence is recovered, 
as the variance and mean are proportional in a binomial distribution - and the cr(Z) — L 
plot is simply recovering the autocorrelation between the mean and the variance of the 
dwarf galaxy Z distributions. This tight sequence in the dwarf galaxies is contrasted by the 
behaviour of the star clusters, which scatter broadly over a range of metallicity spreads. The 
star clusters do not show a correlation in (t{Z) — L since there is no metallicity-luminosity 



relation for GCs in the Milky Way (c.f., Bailin & Harris 2009). 



^We consider all stars in the M54 and Sgr region from this study and note that due to 
the extreme difficulty in associating stars simply to M54 or Sgr (see Carretta et al. 2010), 
the values in this table should be considered upper limits 



- 14 - 



Clearly cr(Z) — L is not an ideal metric for differentiating star clusters and DGs, and 
Figure [2] offers a robust alternative in — Z . This metric does not require knowledge of 

the host object's luminosity, or mass, but simply a distribution of metallicities. While useful 
in classifying low luminosity systems, the binomial nature of the metallicity distributions 
and the Z — cr{Z)^ correlations may also allow for information on chemical enrichment 
processes to be obtained when compared to suitable models. 

5. Discussion 

We have demonstrated that due to the binomial nature of the linear metallicity 
distributions for star clusters and dwarf galaxies, separating them in the (t{Z) — L plane 
is difficult since each class of object has a different luminosity-metallicity relationship. A 
much cleaner distinction between star clusters and DGs is found when plotting the intrinsic 
spread in linear metallicity (characterized by the intrinsic variance) as a function of its 
mean metallicity, cr(Z)^ — Z. In this parameter space, dwarf galaxies form a tight sequence 
over nearly 7 orders of magnitude in luminosity, and separate out cleanly from both open 
clusters and globular clusters. Furthermore, the correlation between Z — <y{ZY is natural 
- a physical result of the sequential buildup to present day metallicities through discrete 
enrichment events in a galaxy that are small in magnitude relative to the total metallicity 
evolution. Understanding the relative scatter and offset in these two sequences offers new 
information on the different modes of chemical enrichment expected to dominate the two 
classes of objects. 
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5.1. Linking Binomial Parameters to Physical Processes in SF 

Having shown that the linear metal distributions are described by binomial 
distributions, it is important to associate the binomial distribution parameters to physical 
processes in the chemical enrichment of a galaxy. As discussed in §3.1, the number of 
successes n and probability of each success q which characterize a binomial distribution 
can be thought of as the number of star formation or enrichment "successes" in a galaxy 
or star cluster, which occur with a characteristic probability. Given the clean separation 
when looking at the dispersion indices, it would be useful to apply analytic chemical 
evolution models to the data that incorporate the binomial form of the Z distributions, 
while providing a link between the statistical properties of the distributions and physical 
processes related to enrichment. There are literature studies that have described galaxy 
chemical evolution using stochastic and inhomogeneous chemical enrichment models (e.g., 
van den Hoek fc de Jong||1997[ |Argast et al.||2000| |Karlsson]|2005[ [Salvadori Ferrara 



2009| ) that can be compared to discrete statistical distributions. In particular, Oey 



(2000) described a modified simple one-zone chemical evolution model based on binomial 
parameters which physically corresponded to the number of enrichment generations, n, and 
the filling factor of chemical enrichment within the galactic ISM, q. In that picture the 
binomial "successes" are physically realized as the successive enrichment of a particular 
volume of gas, which provides an excellent physical link between chemical evolutionary 
processes and the binomial Z distributions found here. 



5.1.1. Oey (2000) Inhomogeneous Chemical Evolution Model 



In the simple closed system that Oey (2000) considered, the inhomogeneous enrichment 



of a galaxy or star cluster proceeds through sequential pollution of regions produced by 
SNe. For n generations of star formation and subsequent enrichment, the galaxy is polluted 
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by bubbles that have some covering fraction q of the whole galaxy. Within a given bubble 
the yield of metals from SNe is assumed to be constant, and therefore the metallicity within 
a given volume is determined only by the size that the bubble sweeps out - as the region's 
metallicity is the result of diluting the pristine ISM with the SN ejecta. 

For a single generation of star formation, the distribution of bubble sizes (and thus 
metallicity of the regions) is well described by a power law, Ns oc with slope (3 = 2 



(Oey & Clarke 1997). The distribution of metals then follows from this distribution of 



bubble sizes: 



......... 



up to some minimum and maximum metallicity. These limits on the Z value of the regions 
in a given generation can alternatively be thought of as corresponding limits on the region 



size, rmin and r^ax- Oey (2000) bounded the minimum bubble size to be that of observed 
supernovae remnants, ~ 25pc in the Milky Way. We relax this constraint in §5.2 to fully 
describe the binomial self enrichment in terms of the physical properties of a galaxy's ISM. 

At any given location in the galaxy, the instantaneous metallicity depends on the 
number of SF generations n and the covering fraction q of the polluting bubbles. This 
location may be polluted by one metallicity region, or it may be overlapped by a second 



generation, or overlapped by up to j regions. Oey (2000) adopted the hkelihood that any 



location is overlapped by j regions is given by the binomial probability: 



n 



P, = (^Jg^(l-gr-^; 1 < j < n. (5) 
In Figure [5] we show a schematic representation of how the number of generations, amount 



of overlap, and range of bubble sizes can be visualized in this framework (taken after Oey 
2003|. 
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The total metallicity distribution function of a host system after n generations is then: 

^ n n 

Ntotiz) = -J2J2 Du-iP,,kNj{z), (6) 
j=i k=j 

Where D is a factor to account for the consumption of gas in the system (see Section 2 of 



Oey 2000 for more details). With these equations from the Oey (2000) model, it is possible 



to create simulated populations for a given n and q using Monte Carlo methods, which has 
the binomial nature of the Z distributions encoded and directly linked to physical processes 
in chemical evolution. 



5.2. Application of the Oey (2000) Model to the a{Zf - Z Diagram 



The Oey (2000) binomial chemical evolution model describes the physical processes 



that control the placement and evolution of star clusters and dwarf galaxies in the cr(Z)^ — Z 
plane. The specific definition of the model parameters could be refined in more complex 
chemical evolution models (e.g., low q could be analogous to describing the fraction of 
ejected/outflowing metals), however here we consider only a general model for our simple 
applications. In addition to the number of SF generations, n and covering fraction g, the 
final metallicity distribution in the model depends on the lowest metallicity of any region, 
Zmim and the range of metallicities /S.Z = Zmax — Zmin- Thus AZ is analogous to describing 
the range of characteristic region sizes, r^m to r^ax found in a galaxy, such that small AZ 
would require a system to have a limited range of bubble sizes (which may be the case when 
the object is intrinsically small as in the case of a star cluster). 

With a choice of initial limits for these bubble sizes/metallicities, the prescriptions in 
the Oey (2000 ) equations permit one to take an object from a starting point on the o"(Z)^ — Z 



diagram, and chemically evolve it an arbitrary amount. The resultant time-integrated 
MDF of the object can be computed through Monte Carlo simulations of mock populations 
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using Equations 4-6. We have computed a grid of simulated MDFs that chemically evolve 
a host from various starting points to quantify the behaviour of objects within this model 
framework. At each end point of a simulated system, the final time-integrated MDF was 
output, and the intrinsic variance and mean of the MDF measured in order to compare 
to the observational data. The 4-dimensional grid of models has a variance and mean 
of the time integrated MDF computed at every point, with the grid points as follows: 
n = [1,2,5,10,25,55] generations; covering fractions of g = [0.1,0.25,0.5,0.75,0.9,1.0]; 
Zmin = [10-^-^, 10-3-0, 10-2-3], and AZ = [0.3, 0.8, 1.6] dex. We present subsets of the model 
tracks in Figure [6] with two parameters fixed, and two varying in each of the panels. 



As noted by Oey (2000), the model is partially degenerate between n and q such that 
the product of those two parameters describes the degree to which a system is chemically 
evolved. This degeneracy can be seen in the top two panels of Figure [6] when comparing the 
metallicity spreads and means. In both cases as the product nq increases, the final object 
moves to a region of higher mean metallicity and metallicity spread, with the binomial 
model doing an excellent job of naturally reproducing the slope of the DG sequence. The 
lower left panel shows that the effect of changing Zmm predictably shifts any starting point 
along the x-axis. The lower right panel illustrates that increasing AZ moves the start point 
to higher metallicity spreads. Notice that for the DGs to have started on their sequence, 
they require a larger range of metallicities within their enrichment regions (which may be 
determined by low mixing efficiencies), or equivalently would need a range of bubble sizes 
larger than what is needed for starting on the GC sequence. This is consistent with simple 



expectation for dwarf galaxies, where intrinsic sizes are larger (> 25 pc; Gilmore et al. 2007 



Tolstoy et al. 2009 ?) than the minimum bubble size corresponding to the smallest region 



adopted by Oey (2000) (<25 pc). 



The variation in starting points as a function of AZ is an important point to consider. 
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as it explains how dirrs with ~ 12 Gyr of continual SF (such as the LMC), can lie on the 
same sequence as dSphs like Segue 1 or Tucana which had all of their stars form in < 1 Gyr 



(Tolstoy et al. 2009). If the number of enrichment events (or equivalently SF duration) was 
the only parameter determining a galaxy's positioning on Figure [6] it would be difficult to 
reconcile how such diverse galaxies could lie on a similar sequence. However, the lower right 
panel illustrates that for specific ranges of AZ a galaxy could start and end on the present 
day DG sequence even with a brief SFH. It should be kept in mind that the application of 



the Oey (2000) model is most useful in a differential sense - i.e., the absolute value of AZ, 
and therefore the physical range of initial DG sizes, will depend somewhat on the yield 
assumed from SNe. 

The model tracks in Figure [6] indicate that objects may never move horizontally - 
and are restricted to never evolve fiatter than the slope of the current DG sequence. The 
restricted changes in the metallicity spread of a system places a fundamental limit on the 
homogeneity of the chemical enrichment in dwarf galaxies. The predictions from the model 
match the observed DGs, and are consistent with the canonical theory of self-enrichment 
dominating in dwarf galaxies. The fact that this holds over 8 orders of magnitude in 
dynamical mass suggests that while the enrichment details differ, the dwarf galaxies 



themselves undergo similar inhomogeneous chemical evolution (or merging; c.f. Salvadori & 



Ferrara 2009 ) processes throughout their lifetime - with their end states well described by a 



simple binomial model. 



5.3. Star Cluster Pre-Enrichment and Dwarf Galaxy Self-Enrichment in the 

Z — o(Z)^ Diagram 



If the chemical evolution in the dwarf galaxies was dominated by self-enrichment, 
whereas the star clusters are pre-enriched, then it is possible to describe both the split 
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in the sequences in Figure |2] and the larger scatter in the star cluster sequence in the 
framework of the inhomogeneous mixing model. In this picture the proto-GCs would have 
had their ISM enriched by the larger environment (a galaxy host) out of which they formed, 
with subsequent SF being too brief, too energetic, or the gas content simply too small 
to allow for continued enrichment within the cluster. The smaller values for at all 

metallicities for the star clusters would be expected then, as they would form stochastically 
from portions of a galaxy's ISM that has been inhomogeneously enriched. On the other 
hand, dwarf galaxies would start their chemical evolution from a relatively un-enriched 
ISM and can undergo a larger change in mean metallicity (increasing (t(Z')^) due to more 
efficient retention of metals and gas in their larger potential wells. An increase in cr{Zy 
for the dwarf galaxies would also be produced if they are fed with un-enriched gas via cold 



flows (Keres et al. 2005 Brooks et al. 2009). It is important to note that the duration of 
star formation is likely not the driving parameter between the two sequences as the SF 
epoch in the clusters is thought to be nearly instantaneous, but the SFHs among the dwarf 
galaxies range from less than a Gyr to a Hubble time and yet they still lie on the same 
sequence. This is most apparent in considering the location of Willman I, one of the least 
luminous systems, which occupies a region very near to the SMC and Fornax. 

The slope of the DG — Z sequence is naturally reproduced as the product of nq 



increases in the inhomogeneous model of Oey (2000) (as discussed in §5.2). The limiting 



slope produced by the models requires either that DGs have been born on the sequence 
with a larger AZ value, or that they start below it but self enrich significantly (panel d. 
Figure |6]). This can set complementary constraints on the initial sizes of DGs early in their 
history, or on the number of stellar generations. If all DGs started at a value of AZ = 0.8 as 
shown in panel d of Figure |6l this may be informing us about the characteristic clustering 



size of the central regions of their DM halos. As discussed by Gilmore et al. (2007), if the 



baryons in DGs are initially distributed such that their characteristic scale length is set 
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by the minimum dark matter clustering length or core size, then this may set q and AZ. 
Alternatively the mid-plane ISM pressure has been suggested to have a dependence on SFR 



(Shetty & Ostriker 2012), which may result in mixing efficiencies which are common to 
most galaxies. In both scenarios chemical evolution proceeds in a singular direction due to 
similar initial conditions or ISM physics, and results in the slope and narrow scatter seen in 
Figure |2] Thus, the small scatter suggests that dwarf galaxies have been dominated by self 
enrichment, consistent with the expectations for their chemical evolution. 

By contrast, self-enrichment among the star clusters would be inconsistent with the 
large scatter in their (t{ZY — Z sequence. Panel c of Figure |6] illustrates the effect of a 
second generation of SF on the GCs. The increased spread in metallicity is far less than 
the scatter shown by the cluster points, suggesting that even with a second generation 
self-enrichment cannot fully explain the star clusters sequence. In addition the direction 
of movement after one generation at such chemically unevolved states is nearly vertically. 
This importance of AZ can also be appreciated when asking if GCs can move diagonally 
along their sequence. The only way to start an object on the GC sequence is to change AZ 
to low values. However this not only starts the evolutionary model at a lower metallicity 
spread, but simultaneously restricts the subsequent generations to increase in dispersion 
faster than mean metallicity. 

A further elegant constraint on the amount of self-enrichment in GCs was demonstrated 



by Bailin & Harris (2009), who showed that the relative spread in cluster metallicity could 
be expressed independent of star formation or metal retention efficiency and only as a 
function of cluster mass: 



-1/2 

This maximum amount of self-enrichment [o(Z^jZ^ that a globular cluster of a given 
mass could undergo provides a strict upper limit to the cluster's ability to enrich, and can 
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be easily compared to the observed ratio of a[Z)/Z. In Figure [t] we show the relative 
metallicity spread a{Z)/Z for the globular cluster points as a function of cluster mass, 
where masses have been computed assuming the luminosity is converted with a constant 
mass-to-light ratio of 2 (i.e., ?). The solid curve shows the maximum metallicity evolution 
{a{Z)/Z) produced via self-enrichment, as a function of star cluster mass. Except for five 
of the clusters, all GCs in the sample show relative metallicity spreads that are too large to 
be produced by self enrichment. This provides an independent check on the ability of GCs 



to self enrich, confirming the results from the application of the Oey (2000) model to the 
offset sequences in Figure [6| 

The inability of self-enrichment to dictate a star cluster's evolution is even more 
striking when one accounts for the fact that most if not all clusters have had two generations 



of SFn Conroy (2011) showed that ejecta from the primordial generation of stars could 



only be produced and retained if the globular clusters were many times more massive than 



their present day stellar mass. Using a diffusive model, Conroy (2011) computed the mass 



enhancement factor {ft = NxMpresent) for the clusters based solely on the strength of their 
observed Na-0 anticorrelations. We have colour coded the globular clusters in Figure [7] 



with the mass enhancement values from Conroy (2011), and it is clear that four of the five 



clusters in Figure [7] which currently lie below the solid self-enrichment limit would have 
had original masses that place them above that curve. This further suggests that globular 
clusters are inconsistent with having been strongly self-enriched. 

If globular clusters are unable to efficiently self-enrich to explain their a{Z)/Z values. 



■^A second generation is required to explain the elemental anticorrelations and morpho- 
logical features in the colour magnitude diagrams of MW GCs, where the second generation 



stars are polluted by the ejecta from AGB stars formed in the first generation (c.f., Gratton 
et aL||2012l). 
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the alternative is that their positioning and scatter on the cr(Z)^ — Z plane is set by 
pre-enrichment. This requires a galaxy to self-enrich inhomogeneously, such that the 
location or time of formation of a GC from that galaxy sets the GC mean metallicity and 
intrinsic spread. In this picture of pre-enriched cluster formation, a parcel of gas within a 
self-enriched galaxy could form a GC at any point during the galaxy's chemical evolution 
- however due to the inhomogeneous mixing implicit in the binomial model, the stochastic 
sampling of the host galaxy's ISM would lead to smaller cr(Z)^ values for the GC than for 
the host DG. This process results in the offset of the star cluster sequence in the o"(Z)^ — Z 
diagram. The slope of the cluster sequence is then a remnant signature of the self-enriched 
history of the galaxy's gas parcel out of which the GC formed, as the clusters are unable 
to evolve diagonally as shown in panel d of Figure |6} This picture of GC pre-enrichment 
is also consistent with the lack of strong correlation between deviations from the best fit 
cluster sequence in Figure [2] and current host properties of GCs - as their pre-enrichment 
levels were dictated by the host galaxy. 



5.4. Future Uses for the oiZY — Z Relationship 

Analysis of the star cluster and dwarf galaxies jointly in the — Z diagram in 

the context of binomial self-enrichent and pre-enrichment may be used to constrain the 
early evolutionary history of the dwarf galaxies and star clusters, but also to study how 
the system of Milky Way star clusters were accreted and/or formed in situ. Under the 



assumption that all GCs originated in galaxies that were self-enriched, then the Oey (2000) 
model in conjunction with cr(Z)^ — Z values place constraints on the possible physical sizes 
or chemically evolved states of the progenitor galaxies. This ability to track the stochastic 



buildup of metallicity over time with the Oey (2000) model, allows one to quantitatively 



describe the formation environment for star clusters - something not possible with usual 
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analytic chemical evolution models. This could be powerful when used with studies that 
place simulated dwarf galaxies on the o'(Z)^ — 2 plane. Fine tuning of feedback or gas 
inflow/outflow efficiencies in those simulations to reproduce the dwarf galaxy sequence 
could result in further insights into the physical processes which dictate the slope and 
scatter in that parameter space. 

6. Conclusions 

We have analyzed the metallicity distributions of a sample of open and globular clusters 
together with dwarf galaxies in the Local Group, focusing on the statistical nature of their 
linear metallicity (Z) distributions. Below is a brief summary of the main findings from 
this analysis: 

• The star clusters and dwarf galaxies show strong correlations between their mean 
linear metallicity (Z) and intrinsic spread in linear metallicity (cr(Z)^). These 
correlations are consistent with a binomial distribution, (i.e., where (t{ZY / Z < 1), 
unfike a{[Fe/H]Y which has no correlation with [Fe/H] 

• Plotting a{ZY against Z for all the objects shows two distinct sequences - one for the 

star clusters and another for the dwarf galaxies. The o'(Z)^ — Z parameter space is 
therefore useful at distinguishing dwarf galaxies from star clusters, without additional 
information (such as L,Mdyn)- We find that o'Cen and M54/Sgr fall on the dwarf 
galaxy sequence, confirming expectations that they represent the nuclei of accreted 
dwarf galaxies. 

• By applying a binomial chemical evolution model to objects in the <j{ZY — Z 
parameter space (e.g., Oey 2000), the slope and offsets of the dwarf galaxy and star 
cluster sequences can be reproduced and explained in terms of the number of star 
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formation generations, the covering fraction, and characteristic size of the enrichment 
events. 

Within this model framework, we illustrate how the inhomogeneity of the self- 
enrichment process for galaxies dictates the slope of their sequence in the ct(Z)^ — Z 
plane. The offset of the star cluster sequence reflects their formation by stochastic 
pre-enrichment from a parcel of gas inhomogeneously mixed within their host galaxy. 
The star cluster slope represents a remnant signature of the self-enriched chemical 
history of the galaxies from which they formed. 

This model confirms expectations that star clusters cannot have been strongly 
self-enriched like the dwarf galaxies, and that pre-enrichment must have been the 
dominant mechanism affecting the star clusters' heavy element properties. This model 
also explicitly describes the connection between the two enrichment modes of galaxies 
and star clusters, in a self consistent and physically motivated framework. 



Future applications of the Oey (2000) model may be able to provide constraints (size, 
ISM density, potential well depth) on the physical environment of the high redshift galaxy 
environments that star clusters are sampling. Also, while the model is simplistic, constraints 
on enrichment and SF physics may be studied by comparing the predictions from high 
resolution SPH galaxy simulations in the Z — cr(^)^ parameter space to observations. 
Through such comparisons, variations produced by changing SF or feedback prescriptions 
in the simulations could be directly related to movement in Z — oiZ^' diagram. 
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Fig. 1. — Sample metallicity distributions for four dwarf galaxies. Left panels show raw 
[Fe/H] distribution, while the right panels shows the distribution of linear metallicities, Z. 
We note that the spread in the linear metallicity distribution is correlated with the mean 
linear metallicity. 
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Fig. 2. — Intrinsic variance cr(Z)^ in Z vs. Z, the mean Z for dwarf galaxies and star 
clusters. Dashed lines represent a linear least squares fit to the dwarf galaxies [black) and 
star clusters (magenta) . Arrows indicate upper limits to the intrinsic dispersions. There is a 
clear separation between the dwarf galaxy and star cluster sequences. Right panel shows the 
same data but with jackknife sampling errors (solid lines) overlaid. Open circles show the 
effects on the the intrinsic variance if the errors reported in the literature were under/over- 
estimated by 50 — 150% in 10% increments for each object. 
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Fig. 3. — Intrinsic spread in [Fe/H] versus average metallicity for dwarf galaxies and star 
clusters. This shows that there is no correlation between the mean [Fe/H] metallicity and 
the metallicity spread o"(Fe) in dwarf galaxies, unlike the correlation seen in Figure [2] 
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Fig. 4. — Dispersion in [Fe/H] and Z versus host object luminosity. Dwarf galaxies (black) 
show a linear relation in a{Z) vs L due to the linear metal fraction being binomially dis- 
tributed and the well known luminosity-metallicity relation. The 5 UFDs from the sample of 



Kirby et al. (2008) (black squares) have been added to show the linear anticorrelation exhib- 
ited below luminosities of L = IO^Lq. Red points with black circles are M54/Sgr and uCen, 
which are both likely associated with former dwarf galaxies. The open black circle represents 



recent estimates of the initial luminosity of the Sgr dSph from the study of Niederste-Ostholt 



et al. (2012). Small arrows indicate upper limits. 
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Fig. 5. — Illustration of the Oey (2000 ) binomial chemical evolution model (adapted from Fig. 



1 in Oey (2003)). Three generations of enrichment (n = 3; green, red, and blue) are shown, 
each have a covering fraction of q. The range in bubble sizes, from r^m to r^ax correspond 
to a range in metallicities {Zmin-, Zmax) within each bubble assuming the constant yield of 
metals within a region is diluted in the corresponding volume. The probability that any point 



is enriched by j overlapping regions is given by the binomial probability P = — 5')" 
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N 
N 



N 
N 

to 




0.01 



0.1 



0.01 



0.10 

Z/Zo 



1.00 



Fig. 6. — Change in the variance and mean for different combinations of n number of 
generations, q covering fraction, minimum bubble metalhcity Zmin-, and range of region 



metalhcities AZ taken from the binomial chemical evolution model of Oey (2000). Lines 



show the computed variance and mean of simulated linear metallicity distributions as the 
models are evolved from n = 1 to n = 55. Black, red and blue points show the same dwarf 
galaxy, globular cluster, and open cluster data from Figure [2] The lower panels show that the 
clusters must have been pre-enriched, unlike the dwarf galaxies which demonstrate strong 
self-enrichment in the upper panels. 
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4.04 



5.07 



6.10 



7.14 



8.17 



9.20 



Mass dependent self-enrichment: Eq. 10 Bailin & Harris 2009 

a^/Z, = 0.059(Mec/10'Mo)"''' 




10" 



10^ 



lO'' 



10' 



10° 



10'' 



Fig. 7. — The relative spread in Z as a function of globular cluster mass. The theoretical 
upper limit to globular cluster self enrichment as a function of proto-cluster cloud mass 



(from Bailin & Harris (2009)) is shown as the solid line. Shading corresponds to total cluster 
mass enhancement at formation relative to the present time, as estimated from the Na- 
O anticorrelations by Conroy| ( 2011[ ). Incorporating the mass enhancement multipliers, all 
clusters but NGC 6739 (darkest blue circle) show metal spreads larger than can be accounted 
for by self-enrichment. 
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Table 1. Spectroscopic [Fe/H] Data for the Literature Sample 



Name 



Study 



R 



AA 



[Fe/H] Method^ 



NGC 


104 


NGC 


6218 


NGC 


6656 


NGC 


5024 


NGC 


288 


NGC 


6752 


NGC 


1904 


NGC 


3201 


NGC 


4590 


NGC 


5904 


NGC 


6121 


NGC 


6171 


NGC 


6254 


NGC 


6388 


NGC 


6397 


NGC 


6809 


NGC 


6838 


NGC 


7078 


NGC 


7099 


M 54/Sgr 


NGC 1851 


IC 4499 



Lane et al. 



Lane et al. 



Lane et al. 



Lane et al. 



Lane et al. 



Lane et al. 



(2011) 



(2011) 



(2011) 



(2011) 



(2011) 



(2011) 



Carretta et a 



Carretta et a 



Carretta et a 



Carretta et a^ 



Carretta et a^ 



Carretta et a^ 



Carretta et a 



Carretta et a 



Carretta et a 



Carretta et a^ 



Carretta et a^ 



Carretta et a^ 



Carretta et a 



Carretta et a 



Carretta et a 



Hankey & Cole 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2009) 



2010) 



2011) 



(2011) 



2241 
242 
345 
180 
133 
437 

58 
149 
122 
136 
103 

33 
147 

36 

144 

157 

39 

85 

65 

103 

124 

43 



10000 
10000 
10000 
10000 
10000 
10000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
40000 
10000 



MRS 
MRS 
MRS 
MRS 
MRS 
MRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
CaT 
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Table 1 — Continued 



Name 



Study 



R 



AA 



[Fe/H] Method^ 



Trumpler 5 
NGC 6791 
NGC 6819 
NGC 6939 
NGC 7142 
NGC 7789 
Berkeley 99 
King 2 
NGC 7044 
Berkeley 39 
M 67 

Melotte 66 
NGC 2141 
NGC 2419 
NGC6205 
NGC 5272 
NGC 6341 
NGC 362 
Terzan 5 
Pyxis 
NGC 2808 
Rup 106 



Warren & Cole 



Warren &; Cole 



Warren & Cole 



Warren & Cole 



Warren &: Cole 



Warren & Cole 



Warren & Cole 



Warren &; Cole 



Warren &; Cole 



Warren & Cole 



Cole et al. 



(2004) 



Warren & Cole 



Warren &: Cole 



Cohen et al. 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2009) 



(2011) 



Kraft & Ivans 



Sneden et al. 



Sneden et al. 



(2004) 



(2003) 



(2000) 



Shetrone & Keane 



Origlia & Rich 



Saviane et al. 



Saviane et al. 



Saviane et al. 



(2012) 



(2000) 



(2004) 



(2012) 



(2012) 



15 
23 
9 

13 
6 

20 
9 
7 

10 
10 
7 
13 
14 
7 

27 

23 

33 

12 

6 

5 

17 
9 



3400 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
10000 
34000 
45000 
45000 
20000 
30000 
25000 
3400 
3400 
3400 



CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
HRS 
HRS 
HRS 
HRS 
HRS 
HRS 
CaT 
CaT 
CaT 



Table 1 — Continued 



Name 



Study 



R 



AA 



[Fe/H] Method'^ 



NGC 5824 
Lynga 7 
NGC 6139 
Terzan 3 
NGC 6325 
NGC 6356 
HP 1 

NGC 6380 
NGC 6440 
NGC 6558 
Pal 7 

NGC 6569 
Terzan 7 
NGC 7006 
NGC 6441 
NGC 6528 
NGC 6553 
WLM 
LMC 
SMC 
Fornax 
Sculptor 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



Saviane et al. 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



Leaman et al. 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2012) 



(2009 



2012) 



Cole et al. 


(20 


05 


); 


Cj 


irrera et al. 


(2008b) 


Carrera et al. 


(2( 


D08a 


); 


Parisi et al. 


(2010 



Battaglia et al. 



Tolstoy et al. 



(2006) 



(2004) 



17 

8 

15 

13 

10 

11 



4 
14 

7 
12 
20 

7 

4 

13 
126 
815 
729 
870 
629 



3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400 
3400,7000 
3400 
3400 
6500 
6500 



CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
CaT 
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Table 1 — Continued 



Name 



Study 



N, 



stars 



R 



AX 



[Fe/H] Method^ 



Sextans 
Carina 
Leo I 
Leo II 

Canes Venitci I 

Draco 

Ursa Minor 

Hercules 

uCen 

Willman I 

Ursa Major I 
Tucana 
NGC 6822 
Segue I 
Bootes I 



Battaglia et al. 



Koch et al. 



Kirby et al. 



Kir by et al. 



Kirby et al. 



Kirby et al. 



Kirby et al. 



Aden et al, 



(2011) 



(2006) 



(2010) 



(2010) 



(2010) 



(2010) 



(2010) 



(2011) 



Stanford et al. 



Johnson et al. 



Martin et al. 



Martin et al. 



(2006) 



(2008 



(2007) 



(2007) 



Praternali et al. 



Tolstoy et al. 



Norris et al. 



Lai et al. 



2009) 



(2009) 



(2001) 



(2010); 



(2011) 



Geha et al. 



(2009) 



180 


6500 


CaT 


327 


6500 


CaT 


827 


7000 


MRS 


258 


7000 


MRS 


170 


7000 


MRS 


298 


7000 


MRS 


212 


7000 


MRS 


11 


20000 


HRS 


378 


2000 


MRS 


246 


13000,18000 


HRS 


8 


7000 


CaT 


10 


7000 


CaT 


14 


3400 


CaT 


26 


3400 


CaT 


4 


5000,7000 


MRS 


25 


2000 


MRS 



'^Method which the the original study used to compute [Fe/H] estimates: Spectrum synthesis 
based techniques on medium (MRS) or high resolution spectra (HRS); empirical calibration based 
on the Calcium II Triplet (CaT) features. 
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Table 2. Intrinsic Metallicity Spreads 



Name 


[Fe/H] 


a{[Fe/H]Y 




Z 




NGC 104 


-0, 


.82 


0.023 


1, 


.73E-01 


1.02E-03 


NGC 6218 


-1, 


.15 


0.061 


8, 


.04E-02 


2.84E-04 


NGC 6656 


-1, 


.95 


<0.197 


1, 


.26E-02 


4.09E-06 


NGC 5024 


-1, 


.89 


0.040 


1, 


.44E-02 


9.26E-06 


NGC 288 


-1, 


.25 


<0.195 


6, 


.22E-02 


1.39E-04 


NGC 6752 


-1, 


.45 


<0.145 


3 


.74E-02 


1.57E-04 


NGC 1904 


-1, 


.54 


<0.035 


2, 


.89E-02 


<5.52E-06 


NGC 3201 


-1, 


.49 


0.048 


3, 


.25E-02 


1.32E-05 


NGC 4590 


-2 


.22 


0.076 


6, 


.12E-03 


1.29E-08 


NGC 5904 


-1, 


.35 


<0.024 


4, 


.52E-02 


<6.13E-06 


NGC 6121 


-1, 


.20 


<0.025 


6, 


.33E-02 


<1.30E-05 


NGC 6171 


-1, 


.06 


<0.043 


8, 


.67E-02 


<7.11E-05 


NGC 6254 


-1, 


.56 


0.052 


2, 


.80E-02 


1.17E-05 


NGC 6388 


-0, 


.40 


0.074 


4, 


.04E-01 


4.38E-03 


NGC 6397 


-1, 


.99 


0.003 


1, 


.02E-02 


4.04E-09 


NGC 6809 


-1, 


.97 


0.044 


1, 


.09E-02 


1.20E-06 


NGC 6838 


-0, 


.81 


<0.033 


1, 


.56E-01 


<1.34E-04 


NGC 7078 


-2, 


.34 


0.011 


4, 


.61E-03 


2.59E-08 


NGC 7099 


-2 


.36 


0.005 


4, 


.42E-03 


3.34E-09 


M54/Sgri^ 


-1, 


.29 


0.445 


1, 


.09E-01 


1.18E-02 


NGC 1851 


-1, 


.14 


0.058 


7, 


.27E-02 


9.74E-05 


NGC 6205 


-1, 


.61 


<0.055 


2, 


.45E-02 


<9.26E-06 
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Table 2 — Continued 



Name 


[Fe/H\ 


a{[Fe/H]Y 


Z 


a{Zf^ 


IC 4499 


-1.15 


<0.107 


7.33E-02 


<3.56E-04 


Trumpler 5 


-0.49 


0.124 


3.60E-01 


7.76E-03 


NGC 6791 


0.39 


0.147 


2.68E+00 


5.75E-01 


NGC 6819 


-0.06 


<0.108 


9.07E-01 


6.32E-02 


NGC 6939 


-0.22 


<0.039 


6.04E-01 


<3.18E-03 


NGC 7142 


-0.19 


<0.139 


6.65E-01 


3.84E-02 


NGC 7789 


-0.23 


<0.103 


6.13E-01 


3.40E-02 


Berkeley 99 


-0.57 


0.034 


2.97E-01 


4.42E-03 


King 2 


-0.42 


<0.133 


3.99E-01 


<1.33E-02 


NGC 7044 


-0.15 


0.181 


7.57E-01 


1.51E-02 


Berkeley 39 


-0.38 


<0.066 


4.20E-01 


<4.03E-03 


M 67 


-0.19 


<0.042 


6.54E-01 


<3.92E-03 


Melotte 66 


-0.49 


<0.049 


3.26E-01 


<1.32E-03 


NGC 2141 


-0.27 


<0.085 


5.46E-01 


1.02E-02 


NGC 2419 


-2.06 


<0.096 


8.86E-03 


3.50E-06 


NGC 362 


-1.33 


<0.034 


4.74E-02 


<1.32E-05 


NGC 5272 


-1.58 


<0.069 


2.65E-02 


<1.49E-05 


NGC 6341 


-2.34 


<0.039 


4.58E-03 


<1.49E-07 


Ter 5 


-0.25 


<0.057 


5.73E-01 


<5.20E-03 


WLM 


-1.28 


0.387 


8.26E-02 


1.62E-03 


LMC 


-0.55 


0.367 


3.83E-01 


3.14E-02 


SMC 


-1.07 


0.349 


1.46E-01 


5.65E-03 
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Table 2 — Continued 



Name 


[Fe 


/H] 


a{[Fe/H]Y 




Z 








Fornax 


-1 


.29 


0.446 


8, 


,40E- 


-02 


2, 


.43E-03 


Sculptor 


-1 


.95 


0.355 


1. 


,77E- 


■02 


1, 


.57E-04 


Sextans 


-2 


.29 


0.387 


1. 


,17E- 


■02 


5, 


.74E-05 


Carina 


-1 


.94 


0.346 


1, 


.68E- 


■02 


1 


.05E-04 


Leo I 


-1 


.46 


0.286 


4. 


,33E- 


-02 


3 


.88E-04 


Leo II 


-1 


.70 


0.362 


2. 


,79E- 


■02 


2, 


.22E-04 


Canes Venitici I 


-2 


.01 


0.446 


2. 


.03E- 


■02 


1 


.69E-04 


Bootes I 


-2 


.59 


0.373 


4. 


.OOE- 


■03 


3 


.89E-06 


Segue I 


-2 


.73 


0.626 


4, 


.93E- 


■03 


1 


.15E-05 


Tucana 


-1 


.97 


0.265 


1. 


,34E- 


-02 


3 


.90E-05 


NGC 6822 


-1 


.20 


0.457 


1. 


.18E- 


-01 


3, 


.78E-03 


Ursa Major I 


-2 


.07 


0.293 


1. 


.24E- 


■02 


4, 


.62E-05 


Willman I 


-1 


.40 


0.324 


5. 


.62E- 


■02 


8, 


.llE-04 


Draco 


-1 


.99 


0.363 


1. 


,72E- 


■02 


9, 


.45E-05 


Ursa Minor 


-2 


.13 


0.339 


1. 


,40E- 


-02 


6, 


.12E-05 


Hercules 


-2 


.70 


0.331 


2. 


,81E- 


-03 


2, 


.83E-06 


wCen 


-1 


.51 


0.253 


3. 


,93E- 


■02 


6, 


.95E-04 


Pyxis 


-1 


.09 


0.246 


9. 


.27E- 


■02 


1 


.57E-05 


NGC 2808 


-0 


.87 


0.106 


1. 


,40E- 


■01 


7, 


.61E-05 


Rup 106 


-1 


.46 


0.062 


3, 


,61E- 


-02 


2, 


.52E-05 


NGC 5824 


-1 


.74 


0.115 


1. 


.93E- 


-02 


2, 


.41E-05 


Lynga 7 


-0 


.57 


<0.035 


2. 


,72E- 


■01 


<4.59E-04 
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Table 2 — Continued 



Name [Fe/H] a{[Fe/H]Y Z a{Z) 



2a 



NGC 6139 


-1.32 


0.069 


4.88E- 


02 


5.50E-05 


Ter 3 


-0.81 


0.070 


1.62E- 


01 


4.59E-04 


NGC 6325 


-1.10 


0.012 


8.06E- 


02 


2.75E-05 


NGC 6356 


-0.51 


0.061 


3.15E- 


01 


2.18E-03 


HP 1 


-1.04 


0.084 


9.53E- 


02 


3.10E-04 


NGC 6380 


-0.51 


0.073 


3.13E- 


01 


2.48E-03 


NGC 6440 


-0.39 


0.077 


4.10E- 


01 


4.94E-03 


NGC 6558 


-0.83 


<0.069 


1.48E- 


01 


5.72E-04 


Pal 7 


-0.49 


<0.082 


3.29E- 


01 


<3.60E-03 


NGC 6569 


-0.82 


<0.088 


1.53E- 


01 


1.12E-03 


Terzan 7 


-0.32 


0.089 


4.93E- 


01 


1.20E-02 


NGC 7006 


-1.32 


0.074 


5.04E- 


02 


6.71E-05 


NGC 6441 


-0.53 


<0.080 


3.03E- 


01 


<3.23E-03 


NGC 6528 


-0.44 


<0.076 


3.64E- 


01 


<3.67E-03 


NGC 6553 


-0.28 


<0.082 


5.34E- 


01 


<9.43E-03 



^Column reports intrinsic spreads - observational measure- 
ment errors have been removed following Equation 1 

'^We consider all stars in the M54 and Sgr region from 
this sample and note that due to the extreme difficulty in 



associating stars simply to M54 or Sgr (see Carretta et al 



2010), the values in this table should be considered upper 
limits 



Note. — Values of upper limits (<) reported where intrin- 
sic spread is smaller than the lowest measurement error of 
an individual star from the original study (see §4). 



